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Abstract 

We consider the exclusive flavor changing neutral current processes B — * Mvv (M = n, K, p, K*) 
in the leptophobic Z' model, in which the charged leptons do not couple to the extra Z' boson. We 
find that these exclusive modes are very effective to constrain the leptophobic Z' model. In the 
leptophobic Z' model, additional right-handed neutrinos are introduced and they can contribute 
to the missing energy signal in B — > M + $ decays. Through the explicit calculations, we obtain 
quite stringent bounds on the model parameters, \Uf b \ < 0.29 and |f7j|'| < 0.61, from the already 
existing experimental data. We also briefly discuss an interesting subject of massive right-handed 
neutrinos, which might be connected with the dark matter problem. 
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I. INTRODUCTION 



Flavor changing neutral current (FCNC) processes have gained continuous interests as 
powerful probes of new physics beyond the standard model (SM) as well as stringent tests 
within the SM. Since FCNC processes are only generated at the loop level within the SM, 
they are suppressed by the loop momentum and off-diagonal Cabibbo-Kobayashi-Maskawa 
(CKM) matrix elements. Therefore, if FCNC exists in the tree level or exists with an 
enhancement factor in the coupling without severe suppression, it would be very sensitive 
to new physics effects beyond the SM. 5-factory experiments such as Belle and BaBar have 
given various opportunities to test those FCNC processes. Among the various terms in the 
effective Hamiltonian responsible for the decay of b quark, the FCNCs are represented by 
QCD penguin and electroweak (EW) penguin operators. There are plenty of QCD or EW 
penguin dominant processes in non-leptonic and semi-leptonic B decays. As an example, 
recent experimental data for QCD penguin dominant B — > Ktx decays appear to be very 
interesting: Branching ratios and direct and mixing- induced CP asymmetries could not be 
consistently explained within the current SM frameworks, for example, QCD factorization 
(QCDF) \ perturbative QCD (PQCD) Q, and soft-collinear effective theory (SCET) 
models [J Although it has been claimed in Ref. 0] that all the current data except for 
the mixing- induced CP asymmetry could be accommodated within the PQCD scheme by 
considering the next-leading-order corrections, there is still room for new physics or new 
mechanism beyond the SM, especially, in the EW penguin sector [^, 0, Q| to explain all the 
experimental data simultaneously. 

The EW penguin contributions have been thoroughly studied in the SM and beyond, 
too: While the EW penguin effect with the photon emission in h — > 57 processes has 
been precisely determined, and already put stringent bounds on the relevant new physics 
scenarios, experimental measurements for the Z-mediated EW penguin effect are now under 
way at B factories through B — ► X££ and B — ► Xvv decay processes. Although it is 
experimentally very difficult to measure, B — > Xvv decay is an extremely good channel for 
the study of Z-mediated EW penguin contribution within the SM and beyond 



121. 
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In many new physics scenarios, an extra U(l) gauge boson is usually introduced after 
the high energy symmetry breaking or as a low energy effective theory jl^l . 14, In 
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general, the extra gauge boson, commonly named as Z' boson, has a possibility to induce 
the FCNCs. In the present work, we focus on the -E^-based leptophobic Z 1 model among 
various possibilities with the extra VjX) gauge boson extension such as the string-inspired 

3 r 

10, 



grand unified theories (GUTs) 13L Il7|. dynamical symmetry breaking models jl8|, string 
inspired models [3], extra dimension models j^j], and little higgs models j^]. It is well 
known that the extra neutral gauge boson can be leptophobic by introducing the kinetic 
mixing term in the -E^-based model. Originally, the leptophobic Z 1 model was introduced to 
explain the Rb-R c puzzle at LEP 22] and anomalous high-i?^ jet cross section at CDF 



Although its original motivation has been disappeared, the leptophobic Z' model certainly 
remains as a viable scenario beyond the SM. The main experimental constraints on the (non- 
leptophobic) Z' boson arise from the electroweak precision measurements at LEP or direct 
searches via Drell-Yan processes at Tevatron. For the leptophobic Z' boson, however, we 
have now only one experimental constraint from the dijet experiments at DO which exclude 



the mass range for 365 GeV < Mz> < 615 GeV |24(. For the future experiments, while the 
CERN Large Hadron Collider (LHC) could serve as a probe of the leptophobic Z' boson via 
PP ~^ 7 + ft; the future e~e + linear collider (ILC) may have no role as far as leptophobia is 
concerned. In this regard, we note that B — > M + ft decays, where M is a pseudoscalar or 
a vector boson, will be a good channel for the probe of the leptophobic Z' model. 

In this paper, we investigate B — > Mvv (M = 7T, K, p, K*) decays in the leptophobic Z' 
model as a possible candidate of new physics in the EW penguin sector. First, we briefly 
introduce the leptophobic Z' model based on the E$ grand unified theory (GUT) in Sec. ITT1 
Sec. IIHI deals with the effective Hamiltonian and decay rates in detail. In Sec. HVI we present 



numerical analysis and discuss its implications. Concluding remarks are also in Sec. |IV| 



II. LEPTOPHOBIC Z' MODEL AND FCNC 



First, we briefly review the leptophobic Z' model. In GUT or string-inspired point of 



251 ] . It is natural that a U{1)' 



view, the E 6 model is a very plausible extension of the SM 
gauge group remains as a low energy effective theory after the symmetry breaking of the Eq 
group. We assume that the E e group is broken through the following breaking chain 

E 6 -> 50(10) x U(l) f 

-> SU(5) x U(l) x x U(l)^ 



^SU(2) L xU(l)xU(l)', (1) 
where U(l)' is a linear combination of two additional U(l) gauge groups with 

Q' = Qi/j cos 9 — Q x sin 9, 

where 9 is the familiar E§ mixing angle. 

The most general Lagrangian, which is invariant under the SM gauge group with an 

sin y ~ ~ 

extra U(l)', allows the kinetic mixing term £ m i x in g = — B^ V Z' ^ between the U(l) 

and U(l)' gauge boson fields. This off-diagonal term can be removed by the non-unitarity 
transformation 

B (l = fl r taa^ l ( 2 ) 
h F cosx 

which leads to the possibility of leptophobia of the physical Z' gauge boson with the E 6 
mixing. Once all the couplings are GUT normalized, the interaction Lagrangian of fermion 
fields and Z' gauge boson can be written as 



where the ratio of gauge couplings A = gQ'/gy, and 5 = — tan \l A QhThe general fermion- 
Z 1 couplings depend on two free parameters, tan# and 5, effectively [26]. The Z' boson can 
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be leptophobic when (Q'+yj-SYgM) = for L and e c simultaneously. L and e c are the lepton 
fields in the 27 representation of E 6 , as shown in Table UJ In the conventional embedding, 



the Z' boson can be made leptophobic with 5 = — 1/3 and tan 9 = J 3/5. In other embedded 
schemes, where the quantum numbers between L and H (or e c , z/ c , and S c ) are switched, 
one can find suitable values for 5 and tan 9 which make the Z' boson leptophobic . 

In the leptophbic Z' model, FCNCs can arise through the mixing between the ordinary 
SM fermions and exotic fermions. Since the mixing of the left-handed fermions may lead 
to the relatively large Z-mediated FCNCs, we only allow the mixing of the right-handed 
fermions. That is, the mixing between dn(sR, 6_r) and Hr can induce the FCNCs when 
their Z' charges are different 16J . After integrating out degrees of freedom of heavy exotic 
fermions and gauge bosons, the FCNC Lagrangian for the b — > q(q = s, d) transition can be 
parameterized as 

A COS v\y 
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TABLE I: Charges of fermions contained in the 27 representation of Eq within the conventional 
embedding 
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where all the theoretical uncertainties including the mixing parameters are absorbed into the 
coupling U% . Then, the experimental constraints on may be obtained by considering 
related non-leptonic B decays such as B — > Kir decays for which one may also find a clue to a 
solution of the recently-popular B — ► i^7T puzzle. But, it is not easy to obtain reliable bounds 
on TJ% because the non-leptonic B decays severely depend on unknown strong phases 27 1 
and, furthermore, the EW penguin amplitude is sub-dominant in most non-leptonic decays. 

It is obvious that if new physics effects are mainly through the EW penguin contributions, 
the EW penguin dominant process can give the most powerful constraints, and it is natural 
to consider the processes such as b — > sj, b — > s(d)££, b — > s(d)vp as appropriate probes of the 
relevant new physics. In the most new physics scenarios such as the minimal supersymmetric 
standard model (MSSM), model parameters are connected to many different observables, 
and therefore, in order to get consistent constraints on the model parameters we have to 
consider all the relevant processes simultaneously. As we will see later, however, one can 
consider only b — ► s(d)uu decays for the leptophobic Z' model and almost free from other 
experimental constraints. 
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III. THEORETICAL DETAILS FOR B -» Mvv (M = vr, K, p, K*) DECAYS 



In the SM, the effective Hamiltonian describing b — > gz/z/ (g = d, s) decays is given by 

H cS (b -> gz/sM^SM) = ^ W ? <^ g 7 M (l - 7 5 )^7m(1 - 7>, (5) 

Z7TV Z 

where Gp is the Fermi constant, a is the fine structure constant and V^- are elements of 
the CKM matrix. The Wilson coefficient Cf is dominated by the short-distance dynamics 
associated with top quark exchange 191 , and has the theoretical uncertainty due to the error of 

nnn 

top quark mass, whose explicit form can be found in literatures |10LI28LI29I]. Therefore, aside 
from the CKM matrix elements, for the b — > qvv decays, there remains only well-controllable 
theoretical uncertainty from C" Q . Although it is not easy to measure experimentally due to 
only single meson with the missing energy in the final state, it can be analyzed through the 
reconstruction of another B meson from the T(AS) — > B + B~ event. 

Because the leptophobia means that the U(l)' charge is zero for all the ordinary left- 
handed and right-handed lepton fields within the SM, the extra leptophobic Z' gives no 
additional contribution to B — > M + Ifi decays if the missing energy is solely due to the SM 
neutrinos. However, v c or S c in the 27 representation of Eq, as shown in Table HI may be 
interpreted as the right-handed neutrino and so there can arise additional contributions to 
the missing energy signal in B —>■ M + E decays. 

Here we take v c as a candidate for the right-handed neutrino and assume that other heavy 
exotic fermions and bosons be integrated out. Then, the effective Hamiltonian responsible 
for b — ■> qv^pR decay is written as 

#cfr(& -> qu R v R ) = . ™ U^Q z v ' R q^{\ + l5 )bu lfl (l + 7 >, (6) 

where 



Qt R = -^f-^, (7) 

and k depends on the details of unification. Here we take k = 1 for simplicity. Then, the 
total branching ratios (BRs) for B —> M + E signals are given by 

B(B -> Mvv) = B(B -> Mv SM v SM ) + B(B -> Mv R v R ). (8) 

First we take all the neutrinos massless, then the differential decay rates for B — ► Mvv (M = 
P, V) in the leptophobic Z' model are given by 

— {B ^P vv)- 2%5 \V tb V tq \ |C 10 | 



dT 

dq 2 



a 2 \V tb V* tq \ 2 \C» w \ 2 \M\,, 



Ism + 



Ml 



2 10 7T 5 



ISM + 



M 



* 2 Q z J R \u%\ 2 

u 2 \V tb V* tq ?\CU 2 



2 \ 2 



X 



SX v q 



V -^- 2 V V (q 2 f + -{ (1 + V^f (V + I2r v q 2 ) A\{q 2 ) 2 



- 2 AltfY - 2 X v (l - r v - q 2 ) Re {A\{ q 2 )A v 2 { q 2 y 



(9) 



(,10) 



where A M = 1 + r 2 M + g 4 - 2r M - 2g 2 - 2r M q 2 with r M = M^/M 2 , q 2 = q 2 /M 2 and P{V) 
in the superscript and subscript stands for a pseudoscalar (vector) meson, respectively. 

Main theoretical uncertainties arise from the hadronic transition form factors for the B 
meson decay into the P(V) meson. However, one can also reduce these uncertainties by 
taking ratios of the BRs, related by SU(3) flavor symmetry or isospin symmetry, such as 
B — > 7i(p)uu and B — > 7i(p)£i> 1(3] . Especially, the corrections to the strict isospin symmetry 
due to mass difference between p ± and p° and form factors are expected to be very small. 



IV. NUMERICAL ANALYSIS AND DISCUSSIONS 



In Table ITT1 we show the theoretical estimates of BRs within the SM and their current 
expe— bounds at B factor . The errors of the SM estimates ha Tah. O are 
mainly due to the hadronic transition form factors and the CKM matrix elements: For the 
form factors, we use recent results from the LCSR model with theoretical uncertainties of 
about 10 to 13 % at zero momentum transfer 
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3: 



fl(0) = 0.258 ±0.031, 
V p {0) = 0.323 ±0.030, 
A{(0) = 0.242 ± 0.023, 
A p 2 (0) = 0.221 ±0.023, 



/*(0) = 0.331 ±0.041, 
V K *(0) = 0.411 ±0.033, 
A?*(0) = 0.292 ±0.028, 
Af(0) = 0.259 ±0.027. 



(11) 
(12) 

(13) 

(14) 
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TABLE II: Expected BRs in the SM and experimental bounds (90% C.L.) in units of 10 



1 

mode 


BRs m the SM 


T71 ■ ill 1 

Experimental bounds 


B -> Kvv 

B — * TTVV 

B -> K*vv 

B — ► /WZ/ 


c 01 +1.11 
°- Ji -1.03 

22+ - 27 
11.15j£* 


< 36 [30] 

< 100 [31] 



The details of the form 
adopted PDG results |34j . 



actors are given in Appendix. For the CKM matrix elements, we 



\V th \ ~ 1, \V ts \ = 0.040 ± 0.003, \Vtd\ = (9.4 ± 4.6) x 10" 



(15) 



Please note that b — > d transition modes receive much larger theoretical uncertainties 
from \V t d\, as can be seen in Table ITT1 The BRs of B — > K^vv decays, which are the 
b — > s transition process, are of an order of 10~ 5 to 10~ 6 . On the other hand, the BRs of 
B — > tt(p)i/D arising from the b —>■ d transition are less by about one order. And as expected 
from larger degrees of freedom in final states, BRs of the vector boson production processes 
are about 2 or 3 times larger than the corresponding pseudoscalar processes. 

Recently the Belle and BaBar Collaborations presented the upper limits on BRs of B — > 
Kvv and B — > -kvv decays j^, [y] . For the K meson production process with the missing 
energy the current experimental bound is about 7 times larger than the SM expectation, 
whereas for the ir meson production the upper bound is much higher by an order of 10 3 than 
the theoretical estimate. Therefore, there is plenty of room for new physics effects in these 
decays at present. The BRs of all the modes are expected to be (more precisely) measured 
at B factories soon. 

Now we investigate additional effects from the leptophobic Z' boson. In Figs. ^ and |2l 
we present our predictions for the BRs in the leptophobic Z' model as a function of the 
effective coupling \U^ b \ in the model. The solid and dotted lines denote the SM predictions 
and the current experimental bounds, respectively. In the leptophobic Z' model, we have two 
parameters, the mass of Z' boson and new FCNC coupling, \U^\. Since the DO experiment 
excludes the mass range 365 GeV < Mz> < 615 GeV [24], we take Mz> = 700 GeV, which is 
also consistent with the mass bound of the conventional non-leptophobic Z' model. Although 
we choose a specific mass for the Z' boson, our analysis can be easily translated through the 
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FIG. 1: Branching ratios for (a) — > K^vv and (b) — > tt^uu, where v can be the ordinary 
SM neutrinos or right-handed neutrinos. 
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FIG. 2: Branching ratios for (a) B^ — > K*vv and (b) B^ — > p^vv, where v can be the ordinary 
SM neutrinos or right-handed neutrinos. 

corresponding changes in the effective coupling |t/J,'| for different Z' masses. Considering 
the current experimental upper bounds, one can extract the following constraints for the 
FCNC couplings from Fig. ^ 



|[/J|<0.29, \U*\<0.61, 



(16) 



for B — > Kvv and B — > -nvv decays, respectively. Compared with the inclusive b — > svv 
decay case we find presently the exclusive mode gives more stringent bounds on the 



leptophobic FCNC couplings. 

As already mentioned, although the exclusive B — ► Mvv decay modes are much easier 
at the experimental detection than the inclusive ones, they have inevitable theoretical un- 
certainties from hadronic form factors. For example, our expectation value for B — > Kvv 
in the SM with recent form factors from LCSR model is about 5.31 x 10 -6 , however, it can 



be changed to ~ 3.8 x 10~ 6 by using different form factors [jj. In order to reduce the form 
factor uncertainties, we take ratios for BiB — > Mvv) to B(B — > Mev) for M = tt, p mesons. 
In the limit of isospin symmetry, and requiring 
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(17) 



(18) 



where the form factor dependencies are all cancelled out. For K^*' meson productions, 
however, since there does not exist such a corresponding semi-leptonic decay mode, one 
cannot adopt the same analysis given above. Instead, considering M p « Mk* and the flavor 
SU(3) symmetry, we obtain the following relation: 



B(B= 



K^vv) 
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= 4.5 x 10" 2 ( 
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V \v ub /v ts \ 

On the contrary to — >• K*vv and B ± —>■ p°e ± v, we cannot use the SU(3) flavor symmetry 
between B ± — > K vv and 5 ± — > 7r°e ± i/ transitions because of the possibly large SU(3) 
breaking due to large mass difference between tt and K mesons. 

Up to now we have assumed all the neutrinos are massless. Aside from the SM neutrinos 
of which masses are indeed very small, the extra right-handed neutrinos in the leptophobic 
Z' model can have sizable masses. Effects of the massive neutrinos will be an interesting 
subject by themselves. Here we briefly introduce the mass effects. In Fig. 01 we present 
the differential BRs as a function of the normalized momentum transfer-squared for m UR = 
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FIG. 3: Differential BRs as a function of the normalized momentum transfer square, q = q 2 /M|, 
in units of 1(T 6 , for (a) B ± -> K±vv, (b) B ± -> n^uv, (c) B ± -> K^vv, and (d) B ± -» p ± uv. 
Here the decay rates in the leptophobic Z' model are normalized to be five times larger than those 
in the SM. 



0, 0.2, 0.5, 1 GeV. We normalized the integrated BRs such that the BRs in the leptophobic 
Z' model are five times larger than those in the SM. It is apparent that the differential BRs 
show the specific features according to the corresponding mass of the right-handed neutrino. 
Especially the sharp rise near the threshold point allows to increase the accuracy of the mass 
measurement of the lightest right-handed neutrino. Unfortunately, if its mass is lower than 
a few hundred MeV, it is very hard to find the difference from the massless neutrino case. 
In general, the additional right-handed neutrinos can have Dirac mass terms coupled with 
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the ordinary SM neutrinos as well as Majorana masses |35j|. Since there is no limit for both 
Dirac and Majorana masses in principle, they are free parameters. One interesting scenario 
is a so-called "z/MSM" model which allows three right-handed (light) neutrinos in 
addition to the SM ones. Assuming that Majorana masses of right-handed neutrinos are the 
order of the electroweak scale or below and Yukawa couplings are very small, this model can 
be accommodated with the present neutrino data. Especially, if masses of the right-handed 
neutrinos are in range of 2 < M UR < 5 keV, they can be good candidates for the warm dark 
matters, and they can explain the neutrino oscillation via the the see-saw mechanism and 
baryon asymmetry of the universe as well j^ . 

To conclude, the FCNC processes with the missing energy signal in B decays have been 
studied in the leptophobic Z' model. In general the direct probe of the leptophobic Z' model 
may be very difficult at ILC as well as LHC because of large hadronic backgrounds. In this 
regard the exclusive FCNC processes B — > Mvv are very adequate to measure new physics 
in the leptophobic Z 1 model with the right-handed neutrinos. This model could be quite 
important in the context of possibly large new physics scenario in the EW penguin sector. 
We also showed that ratios of BRs can reduce the large hadronic uncertainty from form 
factors. The differential BRs are very useful if the right-handed neutrinos have the sub-GeV 
masses. We also noted that the right-handed neutrinos could be accommodated with the 
"z/MSM" scenario. 
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TABLE III: The parameters for form factors in LCSR model 



32, 



where m& = 4.8 GeV is fixed. 
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APPENDIX : HADRONIC FORM FACTORS 



The hadronic matrix elements for B — > Pvv (P = it, K) decays can be parameterized in 
terms of the form factors f+{q 2 ) and f^(q 2 )'- 



(P(p 2 )|g 7 M (l-7 5 )&|£(Pi)> = P^f^ + q'fW) 



(20) 



where p = pi + P2 and q — Pi — p 2 - For B —>■ Vvv (V = p, i^*)decays, the hadronic matrix 
element can be written in terms of four form factors as 

(V(p 2 ,e)\qY(l-l 5 )b\B( Pl )) 



(7 W v (q^ 



e p p 2 q 



■(£*■ q)(pi+p 2 ) 



Mb + My 



i\s*(M B + M v )A v 1 (q 2 ) 



M Mb + M\ 



with A^{0) = A%{0) and A^{q 2 ) 



M B + M v _ AX{q2 



-q^.q) 2 -^{Al(q 2 )-AUq 2 ))\ (21) 
AW)- 



M B -M v v 2 



2 My 1 ' ' 2My 

For the numerical analysis, we follow the theoretical estimates from LCSR model 



32 



The transition form factors are parameterized by three independent fit parameters and one 



13 



resonance mass. Form factors used in the numerical calculation are given by 



; + (g2) ~ 1 - aVml + 1 - oVml ' (22) 



f K ( a 2 ) = ^ + T -l m\ 

J+[q 1 l-q 2 /m 2 R ^ (l-g 2 /m|) 2 ' 1 ' 



fW) = (24) 



F ^ = (for .F = V or Aq) , (25) 



-W) = ; "I, 2 , (26) 



1 - g 2 /m4 (1 - q 2 /mi t ) 2 



where values of parameters are shown in Table ITTT1 and /_ = (/q — /_ 



P - ( fP fP\M 2 B -Ml 
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